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ABSTRACT 

We present an analysis of seven flares detected from five RS CVn-type binaries (UZ 
Lib, a Gem, A And, V711 Tau and EI Eri) observed with XMM-Newton observatory. 
The quiescent state X-ray luminosities in the energy band of 0.3-10.0 keV of these stars 
were found to be lo^o ^-ao.g ^^.g g-i^ rpj^g exponential decay time in all the sample 
of flares range from ~ 1 to 8 hrs. The luminosity at peak of the flares in the energy 
band of 0.3-10.0 keV were found to be in the range of 10^°-* — 10^^-* erg s~^. The 
great sensitivity of the XMM-EPIC instruments allowed us to perform time resolved 
spectral analysis during the flares and also in the subsquent quiescent phases. The 
derived metal abundances of coronal plasma were found to vary during the flares 
observed from a Gem, V771 Tau and EI Eri. In these flares elemental abundances 
found to be enhanced by factors of 1.3 — 1.5 to the quiescent states. In most of the 
flares, the peak temperature was found to be more than 100 MK whereas emission 
measure increased by factors of 1.5 - 5.5. Significant sustained heating was present in 
the majority of flares. The loop lengths (L) derived for flaring structure were found to 
be of the order of 10^""^^ cm and are smaller than the stellar radii (R^) i.e. L/R*< 1. 
The flare from a Gem showed a high and variable absorption column density during 
the flare. 

Key words: stars: activity stars: eoronae stars: flare stars:individual: RS CVn 
star:magnctic field - X-rays: stars 



. 1 INTRODUCTION 

, Solar/stellar flares are known to be a manifestation of the 
■ reconnection of magnetic loops, accompanied by particle 
I beams, chromospheric evaporation, rapid bulk flows or mass 
ejection, and heating of plasma conflned in loops (e.g. Den- 
nis & Schwartz 1989; Martens & Kuin 1989). Observation- 
ally, flares reveal themselves across the electromagnetic spec- 
trum, usually characterized by a rapid increase of radiation 
reaching a peak value that may occur at somewhat difl^er- 
ent times in different wavelength bands (e.g. Osten et al. 
2000) , followed by a gradual decay. Solar X-ray flare classifl- 
cation schemes (Pallavicini et al. 1977) distinguish between 
'compact' flares in which a small number of magnetic loops 
lighten up on a time scale of minutes, and 'two- Ribbon' flares 
evolving on a time scale of several hours. The latter class 
involves complex loop arcades anchored in two roughly par- 
allel chromospheric Ha ribbons. These ribbons define the 
footpoints of the loop arcades. Such flares are energized by 
continuous reconnections of initially open magnetic flelds 
above a neutral line at progressively larger heights, so that 
nested magnetic loops lighten up sequentially, and possibly 
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also at different times along the entire arcade. The largest so- 
lar flares are usually of this type. Analogy of solar flare clas- 
siflcation schemes to the stellar case is problematic as most 
stellar flares radiate several orders of magnitude more energy 
than a solar flare (Giidel & Naze 2009). However, a compre- 
hensive EXOSAT survey of stellar flares by Pallavicini et al. 
(1990) revealed examples of two different type of flares: (1) 
impulsive flares which are like compact solar flares, and (2) 
long decay flares which are like two-ribbon solar flares. The 
impulsive flares have total energy release of 10'^'' erg s^^, and 
typically last for less than an hour. These are very likely to 
be conflned to a single loop. The long decay flares are much 
more energetic (10^^ erg s~^), of longer durations (^ 1 h) 
and are believed to release their energy in an arcade of loops. 

The flares produced by stars in RS CVn-like binaries are 
usually detected only in the UV or X-rays. These UV and 
X-ray flares show extreme luminosities and very hot tem- 
peratures (> lOMK). Flares from these stars present many 
analogies with the solar flares. However, they also show sig- 
niflcant differences, such as the amount of energy released. 
Modeling the dynamic behavior of X-ray flares allows us to 
constrain the properties of magnetic loops in ways that can- 
not be done from analysis of quiescent eoronae that provide 
only a spatial and temporal average over some ensemble of 



2 Pandey & Singh 



structures. Therefore, study of stellar flares is a valuable 
tool for understanding stellar coronae, as these are dynam- 
ical events, which therefore embody different information 
than available from static, time-averaged observations (e.g. 
Pandey & Srivastava 2009; Pandey & Singh 2008; Favata et 
al. 2005; Reale et al. 2004). Analysis of flares in cool giants 
and subgiants gives us an opportunity to probe the structure 
of corona, which in principle can be significantly different 
from the dwarf star coronae. In the cool giants and subgiants 
the gravity is considerably lower than for the cool main- 
sequence stars, therefore, yielding a larger scale height, and 
possibly allowing very extended coronae to develop (Ayres 
et al. 2003). 

In a previous paper (Pandey & Singh 2008), we have 
examined the characteristics of X-ray flares in G-K dwarfs. 
The flares in these G-K dwarfs are similar to the solar ar- 
cade flares, which are as energetic as M dwarfs but are much 
smaller than the flares observed in dMe star, giants and pre- 
main-sequence analogous. In this paper, we present an anal- 
ysis of archival data obtained from XMM-Newton observa- 
tions of five active cool giants and subgiants in RS CVn-type 
binaries, namely UZ Lib, a Gem, A And, V711 Tan and EI 
Eri with the aim to extend the understanding of the spectral 
and the temporal characteristics of X-ray flares in them. 

The paper is structured as follows: in section 2, we give 
the basic properties of stars in our sample, in section 3, we 
describe the observational data sets that we have analysed 
and the methods of data reduction, section 4 contains our 
analysis and results and in section 5 we present our discus- 
sion and conclusions. 



2 BASIC PROPERTIES OF THE SAMPLE 

The basic parameters of the sample of RS CVns are given 
in Table [T] and described below. 

2.1 UZ Lib 

UZ Librae (=BD -08° 3999) contains a giant primary in a 
close binary system with a period of 4.76 days. It is among 
the most active binary systems of the RS CVn class (Bopp 
& Stencel 1981). Grewing et al. (1989) presented lUE spec- 
tra to discover an ultraviolet excess which they attributed 
to a companion star with the effective temperature of 8000 
°K and radius of 1 R©. Olah et al. (2002) rediscussed the 
absolute dimensions of the system and concluded that the 
secondary is most-likely a low-mass cool star and that the 
ultraviolet excess can be explained by the active chromo- 
sphere of the primary. 

2.2 CT Gem 

a Gem (=HD 62044) is a long period (19.6045 days) RS 
CVn-type binary containing an active Kl III type, red giant 
primary. The angular diameter of the giant primary is 2.3 
mas with a corresponding stellar radius of 9.3 Rq (Nordgren 
et al. 1999). The secondary in a Gem has not been detected 
at any wavelength, but restrictions to its mass and the low 
luminosity suggest that it is most likely a late-type main se- 
quence star of under one solar mass (Duemmler et al. 1997). 
It is well studied system at all wavelengths (e.g. Singh et 



Table 1. General properties of stars in the sample 



Name 


Spectral Type 


V 


Period 


Distance 






(mag) 


(days) 


(pc) 


UZ Lib 


KOIII 


9.11 


4.75 


140.3 


CF Gem 


Kim 


5.40 


19.604 


37.5 


A And 


G8III-IV 


3.82 


54.07 


25.8 


Vni Tau 


K1IV/G5V 


5.90 


2.84 


29.0 


EI Eri 


G0IV/M4-5V 


6.95 


1.95 


56.2 



al. 1987; Osten & Brown 1999; Padmakar & Pandey 1999; 
Brown & Brown 2006; Nordon et al. 2006.). 

2.3 A And 

A And (=HD 222107) consists of a G8 lE-IV primary with 
a rotation period of 54 days and an unseen companion in 
an orbit with an orbital period of 20.5 days (Strassmeier et 
al. 1993); the system is therefore asynchronous. It is one of 
the brightest long-period RS CVn's at X-ray wavelengths. 
In the absence of any firm clues as to the identity of the 
secondary, the primary component in A And is thought to 
be wholly responsible for the observed X-ray emission. 

2.4 VTll Tau 

V711 Tau (=HD 22468) is one of the most extensively stud- 
ied members of the RS CVn class of magnetically active close 
binaries. It is a detached, non-eclipsing SB2 system, consist- 
ing of a Kl subgiant star, filling about 80% of its Roche 
lobe and a G5 dwarf star on a 2.84 day orbit (Fekel 1983; 
Donati 1999). The G5V component is also rapidly rotat- 
ing, but its convection zone is less extended, which makes it 
much less active than the subgiant. Therefore, the magnetic 
activity signatures of V711 Tau are dominated by the KlIV 
star. V711 Tau has been extensively studied in the optical, 
the ultra-violet and the radio, while its extreme ultra-violet 
(EUV) and X-ray emission was analyzed in the context of 
surveys of coronal emission from RS CVn systems (Majer et 
al. 1986; Pasquini et al. 1989; Schmitt et al. 1990; Griffiths 
& Jordan 1998). Rotational modulation in the RGS X-ray 
light curve of V711 Tau has been detected by Audard et 
al. (2001), with a minimum flux occurring near the phase 
when the G5 star is in front (phase = 0.5), consistent with 
a previously reported correlation between binary phase and 
X-ray intensity by Agrawal & Vaidya (1988). 

2.5 EI Eri 

EI Eridani = HD 26337 is a rapidly rotating (vsini = 
51kms^) active, noneclipsing, single-lined spectroscopic bi- 
nary star. It consists G5IV primary and M4-5V secondary 
with orbital period of 1.947232 days (Washuettl et al. 2009). 

3 OBSERVATIONS AND DATA REDUCTION 

Stars in the sample were observed with the XMM-Newton 
satellite using many different set-ups of the detectors. The 
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Table 2. Log of observations with the XMM-Newton 



Star 


Rev. 


Instrument 


Start time 


Exposure 


offset 


Name 




(mode, filter) 


date (UT) 


Time(s) 


(') 


UZ Lib 


127 


MOSl(PFW,medium) 


2000-08-19(04:53:55) 


22616 


0.148 






MOS2(PFW, medium) 


2000-08-19(11:09:59) 


22563 








PN(PFW, Medium) 


2000-08-19(05:34:59) 


23567 






210 


MOSl(PFW, Medium) 


2001-01-31(19:30:24) 


27336 


0.095 






MOS2(PFW, Medium) 


2001-01-31(19:30:23) 


27347 








PN(PFW,Medium) 


2001-01-31(20:11:21) 


24896 




(T Gem 


243 


MOSl(PFW, thick) 


2001-04-06(16:40:29) 


49225 


0.068 






MOS2(PFW, Thick) 


2001-04-06(16:30:52) 


49920 








PN(PSW, Thick) 


2001-04-06(16:46:36) 


55811 




A And 


208 


MOSl(T,thick) 


2001-01-26(20:20:58) 


27944 


0.061 






MOS2(PPW, Thick) 


2001-01-26(20:10:31) 


28541 








PN(PSW, Thick) 


2001-01-26(20:24:00) 


28200 




V711 Tau 


310 


MOSl(PPW, Thick) 


2001-08-18(03:32:01) 


11911 


0.228 






MOSl(PPW,Medium) 


2001-08-18(07:25:52) 


11682 








MOS2(PPW,Thick) 


2001-08-18(03:32:01) 


11935 








MOS2(PPW,Mcdium) 


2001-08-18(07:25:51) 


11684 








PN(PSW,Mcdium) 


2001-08-18(03:47:57) 


25016 






221 


MOSl(PFW,Thick) 


2001-02-22(16:40:45) 


9997 


0.092 






MOSl(PPW,Thick) 


2001-02-22(19:36:35) 


23828 








MOS2(PFW,Thick) 


2001-02-22(16:40:43) 


10000 








MOS2(PPW,Thick) 


2001-02-22(19:36:33) 


23800 






27* 


PN(PSW,Mcdium) 


2000-02-01(09:17:29) 


14603 


1.936 


EI Eri 


1769 


MOSl(PPW,Medium) 


2009-08-06(21:00:31) 


37613 








MOS2(TIMING) 


2009-08-06(21:00:26) 


37361 








PN(PSW,Mcdium) 


2009-08-06(21:05:40) 


37471 





Note: FPW: FullPrime Window PPW: PrimcPartial Window, SW: Small window, T: timing, PSW: PrimeSmall Window 

* No MOS observation 



XMM-Newton satellite is composed of three co-aligned X- 
ray telescopes (Jansen et al. 2001) which observe a source si- 
multaneously, accumulating photons in three CCD-based in- 
struments: the twin MOSl and M0S2 and the PN (Striider 
et al. 2001; Turner et al. 2001), and with all three detectors 
constituting the European Photon Imaging Camera (EPIC 
). The EPIC instruments provide the imaging with a good 
angular resolution [point spread function (PSF) = 6 arc- 
sec full width at half-maximum (FWHM)] and spectroscopy 
in the energy range from 0.15 to 15 keV with a moder- 
ate spectral resolution {E/AE « 20 — 50). Exposure time 
for each star was in the range of 10 -56 ks. A log of ob- 
servations is provided in Table [51 Star V711 Tau was ob- 
served by XMM-Newton in 31 occasions. In 10 occasions 
the exposure times were less then 4.0 ks, therefore, we have 
not considered these observations for further analysis. In 12 
occasions only RGS observations were made. In one occa- 
sion (obs ID:0116710901) the only MOS CCDs were on and 
data is heavily piled up, therefore, is not useful for further 
analysis. In rest of eight occasions, only in four occasions 
flares were observed (see Table 2 for detail on observations), 
in other four occasions (obs ID: 0116340601, 0116320801, 
0129350201, 0134540601) light curves appear to be constant. 
However, during the obs ID 0129350201 a small part of rise 
phase of a flare was observed at the end of the observations. 
The flare observed during the obs ID 0117890901 was anal- 
ysed by Audard et al. (2001). EI Eri was also observed in 
four occasions. 

The data were reduced with standard XMM-Newton 
Science Analysis System (SAS) software, version 10.0 with 
updated calibration files. The preliminary processing of raw 



EPIC Observation Data Files was done using the epchain 
and EMCHAIN tasks which allow calibration both in energy 
and astrometry of the events registered in each CCD chip 
and combine them in a single data file for MOS and PN de- 
tectors. The background contribution is particularly relevant 
at high energies where coronal sources have very little flux 
and are often undetectable. Therefore, for further analysis, 
we have selected the energy range between 0.3 and 10.0 keV. 
Event list files were extracted for each observations using the 
SAS task EVSELECT. Data from the three cameras were indi- 
vidually checked for the time intervals with high background 
when the total count rates (for single events of energy above 
10 keV) in the instruments exceeded 0.35 and 1.0 counts 
s~^ for the MOS and PN detectors, respectively. No data 
were found to be affected by high background proton fiaring 
events. The epatplot task was used for checking the exis- 
tence of pile- up in the inner regions of the X-ray images of all 
stars. Only the MOS data of the star V711 Tau were found 
to be affected by pile-up. X-ray hght curves and spectra of all 
target stars were generated from on-source counts obtained 
from circular regions with a radius 40-50 arcsec around each 
source. However, for the star V711 Tau, X-ray light curves 
of MOS detector were extracted with events taken from an 
annulus of 42 arcsec with an inner radius of 12 arcsec to 
avoid the pile-up effect in the inner region. The background 
was taken from several source-free regions on the detectors 
at nearly the same offset as the source and surrounding the 
source. 
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Figure 1. MOS and PN Light curves of the sample stars. The MOS light curves are scaled by adding 33, 2, 16 and 3.5 counts s ^for 
flares F2, F3, F6 and F7, respectively. 
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4 ANALYSIS AND RESULTS 
4.1 X-ray light curves 

X-ray light curves of the five stars in the energy band 0.3 
-10.0 keV are shown in Figured] The MOS and PN light 
curves are represented by solid triangles and circles, respec- 
tively. The temporal binning of light curves is 200 s for all 
sources. The light curves of the stars UZ Lib (Fig. Hb), o" 
Gem (Fig. [!};), A And (Fig.[TJi), V711 Tau ^ to g) and EI 
Eri (Fig.[T]h) show variability on a time scale of hrs, most of 
which resembles flaring activity. The customary definition of 
a flare is a signiflcant increase in intensity, after which the 
initial or quiescent level of intensity is reached again. Such 
patterns are observed in the light curves of the sample stars 
shown in Fig. [T] where the "flare regions" are represented 
by arrows and marked by Fi, where i=l,2...7 refers to the 
flare number. No flare like features were observed during 
the satellite revolution of 127 in the light curves of the star 
UZ Lib (Fig. [1^), therefore, assumed to represent its quies- 
cent state. The mean quiescent level count rates of the other 
sources are taken from regions that are free from flares and 
are marked by Q in the Fig. [1] For the stars a Gem and 
A And, the quiescent states were not observed. Using the 
RGS data of flare F2 of a Gem, Nordon et al. (2006) found 
an overall 25 per cent flux increase from its quiescent state, 
which was observed 15 months earlier by Chandra X-ray ob- 
servatory. Therefore, we have assumed that the flux during 
the quiescent state of a Gem was 25 per cent less than that 
of the peak flare flux. The flux during the quiescent state of A 
And was derived from ASCA and ROSAT data by Ortolani 
et al. (1997) and is ~ 31 per cent less than the peak value of 
the flare F3. The duration of the flares and the flux values 
during quiescence and at flare peak are listed in Table [3] for 
all the seven flares observed . The peak and quiescent state 
count rates as given in Table [3] are converted into fluxes by 
usmg the WebPIMM^I, where we have assumed a coronal 
plasma temperature of 1 keV. The flux at the peaks of flares 
was increased by 26 - 56 per cent from the quiescent states. 

Peak count rates (Ao), and 1/e rise (r^) and decay (rd) 
times have been determined from the least-squares flt of the 
exponential function c{t) — Aoe'*"^''''^'"''' + q to the count 
rate histogram. Here, to is the time at flare peaks and q is 
count rates during the quiescent state. The fitted values are 
given in Table [3] The exponential decay times were found 
in the range of 0.8 to 9 hrs. However, the 1/e rise times of 
the observed flares were found in the range of 0.3 to 6 hrs. 
Flare F2 of a Gem was found to be the strongest flare, where 
total energy (Etot) released was 4.24 x 10^^ ergs. All other 
flares were weaker by at least a factor of 10 (see Table [3)). 
The MOS data have less count rates than the PN but have 
similar characteristics to that of PN; therefore, for further 
analysis in this paper we use only the PN data. However, 
for the flare F5 from V711 Tau only MOS data were used 
as PN observations were piled up. 

X-ray light curves were also divided separately into 
three energy bands namely soft (0.3-0.8 keV), medium (0.8- 
1.6 kev) and hard (1.6-10.0 keV). The boundaries of the se- 
lected energy bands are chosen as the line- free regions of the 
low-resolution PN spectra. Hardness ratios HRl and HR2 
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as deflned by the ratio of medium to soft band, and hard 
to medium band count rate, respectively, were also derived. 
The soft, medium and hard bands intensity curves, and the 
hardness ratio (HRl and HR2) curves are shown in the sub- 
panels running from top to bottom in Fig.[2]for all the stars. 
The light curves in individual bands show significant vari- 
ability in all sources. The peak fiare fiux, quiescent state 
flux, and 1/e rise and decay times in the three energy bands 
are given in Table g] For the flares Fl, F5, F6 and F7 the 
1 /e decay times in the soft band were found to be more than 
the decay times in the medium and the hard band. However, 
this difference between 1/e decay times in the different en- 
ergy bands is within 1.6cr level. The 1/e decay time for the 
flare F2 was found to be more in the hard band than in the 
soft and the medium bands. The light curves in different en- 
ergy bands, for a Gem, show different behavior. In the soft 
band the light curve shows a flat top, in the medium band 
it shows a narrower top, whereas, in the hard band the light 
curve of a Gem hardly shows the rise phase of the flare. 

The ratio of flare peak flux to the quiescent state flux in 
the hard band was found to lie between 1.4 and 2.5, which is 
more than that of the similar ratios in the medium and the 
soft band (see Table [4]) for all flares. In terms of the peak 
flux, the flare F2 of a Gem was stronger in the medium band 
than in the soft and hard bands. The variation in the hard- 
ness ratio is indicative of changes in coronal temperature. 
An X-ray flare is well deflned by a rise in the temperature 
and subsequent decay. Excluding flare F2, HRl was found 
to be constant during all the flares . However, for most of 
the flares HR2 varied in the similar fashion to their light 
curves. This implies an increase in the temperature at flare 
peak and subsequent cooling. In the case of flare F2, the HR 
variation follows the trend of the hard band light curve. No 
variation was observed in the HR2 of the flare detected from 
UZ Lib. 

4.2 Colour-colour diagrams 

Hardness ratios plotted in the form of colour-colour (CC) 
diagrams can reveal information about spectral variations 
and serve as a guide for a more detailed spectral analysis. 
For this purpose, the light curves and the hardness ratio 
curves of the flares were divided into different time segments 
covering their rising and decaying phases as shown in Fig. 
[T] where Ri(i=l,2..) represents the rising phase, Di(i=l,2,..) 
represents the decay phase. The HRl and HR2 values 
obtained from the PN data for each of the flare segments 
are plotted in Figs. O (a-f). To understand the observed 
behavior of HRl and HR2 in terms of simple spectral mod- 
els, we generated the soft, medium and hard bands count 
rate using the "fakeit" provision in the XSPEC, using the 
most recent Canned Response Matrices downloaded from 
http: / / xmm.vilspa.esa.es / external /xmm_S W_cal / calib /epic] 
_files.html, and simple two temperature plasma models. 
The two temperature APEC models were mixed such that 
the emission measure were in the ratio EM1/EM2 ranges 
from 0.5 to 2.0. The abundances (0.1 to 0.3) and hydrogen 
column density (10^° cm~^) were kept fixed for each set of 
the EM1/EM2. Further by keeping kTi fixed at different 
values ranging from 0.35 to 0.9 keV in steps of 0.05 and 
varying kT2 from 1.0 to 5.5 keV such that kT2 >kTi, we 
generated count rates in the soft, medium and hard bands. 
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Figure 2. PN light curves at three bands soft (0.3-0.8 keV), medium (0.8-1.6 keV), and hard(l. 6-10.0 keV) and hardness ratio HRl and 
HR2 curve, where HRl=medium/soft and HR,2=hard/medium during flares from program stars 



These predicted count rates were used to generate the 
hardness ratios. The famihes of curve thus generated are 
over-plotted on the data points in Fig. [3] to understand the 
CC diagram where kT2 increases from bottom to top of 
each curve shown by arrows. 

For flare Fl, most of the flare segments intersect the CC 
diagram near kTl= 0.6 and 0.7 keV, indicating that the tem- 
perature during the flare remains constant. The best match 
CC diagram for the flare Fl is shown in Fig. [3] (a), where we 
assumed EM1/EM2 = 0.5 and Z = 0.2^©. Fig. [3(b) shows 
the CC diagram for a set of EKh/EKh = 1 and Z = 0.3^0 



for the star a Gem, where the flare peak part D2/D3 inter- 
sects the generated CC curves near kTl = 0.8 keV, however, 
the bottom part Dll intersects the curve near kTl = 0.75 
keV. This imphes that at the flare top plasma temperature 
is high than the flare bottom. A similar pattern was also 
seen for all the flares observed (see Figs. [3j:-f), where the 
top of the flare is located near the CC curve for the high 
temperature component whereas the colours near the end of 
the flare intersects the CC curve corresponding to the lower 
temperature. It appears that the quiescent states of the stars 
are cooler than the corresponding flare states. The temper- 
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Figure 3. Colour-colour diagram for the various time segments as shown in Fig. [T] Here Q represents quiescent state, Ri(i = 1,2,3) 

represents rise phase and Di(i=l,2, ) decay phase of the flare. The curves overplotted are for model simulation where two-temperature 

plasma is assumed to predict the colours. The increasing directions of kT2 is shown by arrows. 



ature is high at the flare peak, and the plasma temperature 
decreases as flare decays. For flare F3, F6 and F7 the simu- 
lated CC diagrams are well matched with the observed one 
for the parameters kTi= 0.4-0.8 keV, kT2 = 1.0 - 3.0 keV, 
Z= O.IZq and EM1/EM2 = 1.0. However, for the the flare 
F5 the observed CC diagram are well matched with the sim- 
ulated families of CC curve for a set of EM\/ EM2 = 0.5 and 
Z = 0.2Zq. 

4.3 The Quiescent coronae 

X-ray spectra for each star during the quiescent state were 
analysed using XSPEC, version 12.5 (Arnaud 1996). Spec- 
tral analysis of EPIC data was performed in the energy band 
between 0.3 and 10.0 keV. Individual spectra were binned 
so as to have a minimum of 20 counts per energy bin. The 
EPIC spectra of the stars were fitted with a single (IT), two 
(2T) and three (3T) temperature collisional plasma models 
known as APEC (Smith et al. 2001), with variable elemen- 
tal abundances. Abundances {Z) for all the elements in the 
APEC were varied together. The interstellar hydrogen col- 
umn density [Nh ) was left free to vary. For all the stars, no 
IT or 2T plasma models with solar photospheric (Anders & 
Grevesse 1989) abundances {Zq) could fit the data, as un- 
acceptably large values of = (where v is degrees of 

freedom) were obtained. For the quiescent state of UZ Lib, 
acceptable APEC 2T flts (x^=1.21) were achieved only when 
the abundances were allowed to depart from the solar values. 
Gondoin (2003) had also analysed the same data of UZ Lib 
and achieved the best fit by using the 4T apec model with 
X? of 1.18. It appears that there is not much improvement in 
while fitting with 4T APEC model. Therefore, we assume 
that quiescent corona of UZ Lib is best represented by 2T 
plasma. For the quiescent state of V711 Tau and EI Eri, only 
3T plasma models with subsolar abundances were found to 
be acceptable. The best-fitting plasma models along with 



the significance of the residuals in terms of ratio are shown 
in Fig. g] for the stars UZ Lib, V711 Tau and EI Eri. Ta- 
ble [S] summarizes the best-fitting values obtained for the 
various parameters along with the minimum smd the 90 
per cent confidence error bars estimated from the minimum 
+ 2. 71 (Arnaud 1996). 
The quiescent states of a Gem and A And were not 
observed during the XMM-Newton observations. The 3T 
plasma models were fitted to the X-ray spectra of each fiare 
segment (for details see section 4.4.2 and 4.4.3) of these two 
stars. In the spectral fitting, all the parameters were left 
free to vary. Figs. [5] and [6] show the variation of tempera- 
tures (Ti and T2) and emission measures (EMi and EM2) 
corresponding to the cool components during the flares for 
the stars a Gem and A And, respectively. The temperatures 
and emission measures corresponding to the cool compo- 
nents are found to be constant during the flare state of both 
stars. Therefore, it appears that the temperatures and emis- 
sion measures corresponding to the cool components are re- 
lated to the quiescent coronae of these stars. The average 
values of the Nh, Z, kTi, kT2, EMi and EM2 for these two 
stars are also given in Table [5] The average value of EM\^2 
and fcri,2 are determined using the following formulae 

EM = - V EMi and kT = V EMikTi 

n ^ EM ^ 

Here, n is total number of flare segments. The derived values 
of quiescent state temperatures and emission measures for 
the star a Gem and A And are found to be similar to that 
derived from the ROSAT data (Shi et al. 1998; Ortolani et 
al. 1997). 

The derived values of Nh from the spectral analysis 
were found to be lower than that of the total galactic Hi 
column density (Dickey & Lockman 1990) towards the di- 
rection of that star. The cool and hot temperatures of qui- 
escent coronae of these stars were found in the range of 0.3- 



Table 4 Parameters obtained from the fitting of an exponential function to the flare light curves observed with PN/MOS detector in the soft, medium and the hard bands. 



Object FN Soft (0.3-0.8 keV) Medium (0.8-1.6 keV) Hard (1.6-10.0 keV) 

Name Td{s) Tr{s) Wp Wq rd(s) r,.(s) F\, F\, Td{s) Tr{s) Fj, 



UZ Lib 


Fl 


29763 ± 2726 




1.3 


0.67 


25038 ± 2365 




1.31 


0.74 


27095 ± 1662 




0.73 


0.29 


a Gem 


F2 


15572 ± 2736 


4557 ± 671 


25.4 


22.5 


17763 ± 2874 


7362 ± 915 


35.0 


29.3 


25167 ± 1068 




25.8 


16.1 


V711 Tau 


F5 


9875 ± 744 


1645 ± 288 


9.5 


6.5 


8921 ± 430 


1433 ± 148 


12.0 


7.3 


7673 ± 336 


1085 ± 77 


10.6 


4.2 




F6 


8731 ± 1737 


6010 ± 879 


16.2 


12.6 


5836 ± 801 


4446 ± 677 


14.4 


11.2 


4323 ± 453 


4279 ± 544 


4.1 


2.6 


EI Eri 


F7 


4201 ± 361 


1227 ± 127 


3.7 


2.7 


3486 ± 247 


1062 ± 76 


9.0 


6.0 


3385 ± 199 


760 ± 55 


5.0 


3.7 



Note: FNMisce is flare Name, Fp and Fq are the fluxes at flare peak and during the quiescent state in units of 10 erg s ^cm ^, respectively, and are flare decay 
and rise times, respectively. 
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Table 3 Flare duration and other parameters obtained from least square fitting of exponential function to the count rate histograms obtained with MOS and PN detectors 
in the energy band of 0.3-10.0 keV 



Object 
Name 



FN 


Flare 








duration 








(ks) 






Fl 


27.0 


0.82 


0.45 


F2 


53.0 


20.97 


15.00 


F3 


> 28.2 


4.54 


3.12 


F4 


3.8 


6.5 


4.8 


F5 


16.6 


11.0 


4.8 


F6 


20.0 


8.67 


4.8 




n n 


Q 71 


1 OQ 



MOS 



PN 



A„ 



Tr 

(s) 



(s) 



A„ 



Tr 



Td 

(s) 



1.80 


133.5 


36.0 


62.3 


800.3 


424. 


12.6 


26.2 


7.4 


26.1 


12.6 


0.4' 




74.8 


12..^ 


26.1 


91.3 


18.2 


1 1 v 


70 fi 


7 1 ( 



UZ Lib 
(7 Gem 
A And 
V711 Tau 



2.08 ± 0.02 
41.2 ±0.12 
13.2 ±0.08 

11.87 ±0.09 
12.4 ±0.1 

11 o -I- n Q 



21802 ± 1059 1.20 



22036 ± 1658 
19587 ± 1034 

1600 ± 127 
6301 ± 521 

1 O^Q -L- QR 



30466 ± 505 



8667 ± 289 
12388 ± 631 

QTKR -I- OOQ 



30.0 
9.2 

7.73 
8.7 

R nc; 



3.28 ±0.04 

84.2 ±0.2 

18.3 ±0.1 
30. ±0.1 

35.5 ±0.1 

on Q -I- n Q 



22012 ± 1571 
19728 ± 995 
1100 ±224 

5190 ± 358 

1 1 OQ -I- '7/1 



32492 ± 523 
2850 ± 521 
6976 ± 169 

QVSQ 4- 1 R? 
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Figure 4. Quiescent state spectra of MOS and PN data with 
APEC 2T/3T model in upper subpancls of each graph, con- 
tributions in terms of ratio arc given in the lower subpanel of the 
each graph. 



0.8 keV and 0.9 - 2.0 keV, respectively. The corresponding 
EM2/EM1 during the quiescent states of these stars were 
found in the range of 1 - 4. 



4.4 Spectral evolution of X-ray flares 

The spectra of the different time intervals shown in Fig. [1] 
were analysed to trace the spectral changes during the flare. 
The temporal evolution of the spectral parameters of the 
flares detected from the targeted stars are shown in Figs. [7] 
- inl and are described below. 



4.4.1 Flare Fl from UZ Lib 

The flare Fl was divided into nine time bins. To study the 
flare emission only, we have performed 3T spectral fits of 
the data, with the quiescent emission taken into account by 
including its best-fitting 2T model as a frozen background 
contribution. This is equivalent to considering the flare emis- 
sion subtracted of the quiescent level, and allows us to de- 
rive one 'effective' temperature and one EM of the flaring 
plasma. Initially, the Nh and Z were left free to vary in 
the spectral fitting of all flare segments and were found to 
be similar to that of quiescent state. Therefore, in the next 
stage of spectral fitting both were kept fix. Fig. [7] shows the 
temporal evolution of temperature and emission measures 
along with X-ray light curves of the flare Fl. The tempera- 
ture was found to be constant during the flare Fl. However, 
the corresponding emission measure was found to decreased 
by a factor of ~ 2.5 along to the decay path. No flrm con- 
clusion can be drawn for the evolution of temperature and 
emission measure as rise and peak phases of the flare were 
not observed. During the flare Fl, the star was found ~ 2 
times more X-ray luminous than its quiescent state (see Ta- 
bles [5] and [6]) . 

4.4.2 Flare F2 from a Gem 

As mentioned in i |4.3l no quiescent state was observed for 
the star a Gem. We have made fourteen time bins for flare 
F2 to study the flare evolution. The temperature and the 
corresponding emission measure of the hottest component 
from the best fit 3T plasma model were found to vary along 
the fiare, therefore, appear to be related to the fiare emis- 
sion. As shown in Fig. [8l^b), the A'^^r increased during the rise 
phase and decreased during the decay phase of the fiare F2. 
The value of Nh reached 2.7 x 10^" cm~^, which is a factor 
of ~ 2 more than the minimum observed value. The temper- 
ature and the Z appear to be highest at the very beginning 
of the flare (Rl) and decrease thereafter (see Figs. [H]c and 
d). The highest value of Z were found to be 0.32Zq, which 
is ~ 1.5 more than the minimum value observed. During 
the flare F2, the temperature evolved before the emission 
measure did. It peaked during the rise phase Rl, whereas 
the emission measure peaked during the decay phase of the 
flare F2. The emission measure was found to vary along 
the flare intensity. The temperature and emission measure 
were increased by factors of ~ 1.8 and ~ 1.3 times, respec- 
tively. The highest luminosity during the flare F2 was 10^^'*^ 
erg s~^, which is three to six times more than that observed 
from EXOSAT (Singh et al. 1987) and ROSAT (Yi et al. 
1997) observations. 

4.4.3 Flare F3 from A And 

In the flare F3, only the rise and the peak phases were ob- 
served. Therefore, an approach similar to that applied for 
flare F2 was applied to the spectral analysis for five time 
bins of this flare. As shown in Fig. [5] (b) and (c), the Nh 
and Z were found to be constant during the rising phase of 
the flare F3. The evolution of the temperature and the emis- 
sion measure corresponding to the hottest component in the 
best flt 3T plasma model are shown in Figs. [9] (d) and (e), 
respectively. The emission measure was found to increase as 
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Table 5. Spectral parameters derived for the quiescent emission of each program star from the analysis of PN and MOS spectra 
accumulated during the time interval 'Q' as shown in Fig. 1, using APEC 2T/3T models. 



Object 




Z 


kTi 


EMi 




EM2 


kT3 


EM3 


Lx 


xl (DOF) 


UZ Lib 
(T Gem* 
A And* 
V711 Tau 
EI Eri 


n q+0.2 
^•■^-0.2 
n Q+0.6 
'^•'-'-0.6 

1-8-0.5 
1 9+0.4 
-^•^-0.4 


'^■^^-0.02 
^•^°-0.02 

^•^°-0.02 

17+0-02 
'^•^'-0.02 


34+0-02 
36+0 01 

42+0 02 
'^■^^-0.02 


o q+0.7 

26.0t«;° 
24.7+3 3 
27.0+3-* 
20.8^2;* 


1 95+0-06 

-, -,0+0.07 
1-18-0.07 
93+0 02 
'^•^'-'-0.02 
95+0 02 
^•^°-0.02 
93+0-02 
'^•^'J-0-02 


S 1+0-3 
94 S~* 
40 5-* ° 
43.815-5 
45.6t|j 


9 c+0.3 
9 c:+0.2 


17.6+2? 


8.69 

5.95 
7.11 


1.21(931) 

1.19(545) 
1.01(655) 



Note: Njj is in unit of 1020 cm 2^ z is the best fit values of abundances, temperatures (kT) are in keV, emission measures (EM) are in 



xl is the minimum reduced and DOF stands for degrees of freedom. 



10 cm ^, and X-ray luminosity (Lx) is in 10 erg s 
* Average values of parameters corresponding to the cool components from the best fit 3-T plasma model to each fiare segment. 
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Figure 5. -Variation of (a) kTi and YT2 and (b)EMi and EM2 
during the flare F2 from a Gem. 
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Figure 6. "Variation of (a) kTi and kT2 and (b)EMi and EM2 
during the flare F3 from A And. 



flare rises. However, the temperature peaked during the flare 
segment 'R2' and decreased as flare rises. It appears that the 
temperature had peaked before the emission measure. The 
luminosity reached a value of 4.6 x 10^" erg s^^, which is 
~ 25% more than that observed ~ 5 years before by ASCA 
and and ~ 10 years before by ROSAT during its quiescent 
state( see Ortolani et al. 1997). 



444 Flares F4, F5 and F6 from ¥711 tau 

Three flares (F4,F5 and F6) were observed during three dif- 
ferent observations of V711 Tau. The flare F4 was a shortest 



duration flare among all observed flares. Therefore, only two 
time bins were made during this flare. However, the flares 
F5 and F6 were divided into 8 and 10 time bins, respec- 
tively. For the flares F4, F5 and F6 of V711 Tau, we have 
fitted the 4T APEC model to the spectral data of each time 
bin, with the quiescent emission taken into account by in- 
cluding its best-fltting 3T model as a frozen background 
contribution. The Nh and Z were left free to vary during 
the spectral fltting. For the flare F4, the values of Nh and Z 
were found similar to that of quiescent state, therefore, were 
not allowed to vary during the spectral fitting. The temper- 
ature and the corresponding emission measure were found 



A study of X-ray flares II. 11 



3.0 




6.0 - 



80.0 



60.0 



S 40.0 




(c) 



5 10 15 20 25 

Time (in ks since .ID 245)941.3127) 

Figure 7. Variation of (a) count rate, (b) temperature and (c) 
emission measure during the XMM-Newton observations of the 
flare Fl from UZ Lib. 



to peaked during the decay phase of the flare F4 (see Table 
[6l. The variation of the temperature, the emission measure, 
Z and Nh along with count rate for the flare F5 and F6 
are shown in Fig. [10] and Fig. 1111 respectively. The A'^^f was 
found to be constant with in la level during the flares F5 
and F6. The values of Nh during flares F5 and F6 were also 
found to be similar to that of quiescent state. The Z were 
found to vary during the flares F5 and F6. During the flare 
F5, the Z increased by factor of 1.5. However, this increment 
is well within 1.7a level from the quiescent state. In case of 
flare F6, Z varied in similar fashion to that of count rates 
and increased by a factor of ~ 1.5 from the quiescent state. 
The emission measure varied along with the rise and decay 
phase of the light curve for the flares F5 and F6 (see Figs. 
llOb .e and lllb .e). However, the temperature peaked before 
the emission measure. The temperature and emission mea- 
sure during the flare F5 increased by factors of ~ 3 and 
~ 2, respectively. However, ratio of temperature and emis- 
sion measure at flare peak and flare bottom for the flare F6 
were found to be 2.8 and 2.6, respectively. The X-ray lumi- 
nosities at the flare peaks were found 1.1, 1.8 and 1.5 times 
more than that of quiescent state for the flares F4, F5 and 
F6, respectively. 
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Time (in ks since ,ID 2452006.1990) 

Figure 8. Variation of (a) count rate, (b) hydrogen column 
density (N//), (c) abundance (Z), (d) temperature (T) and (e) 
emission measure (EM) during the XMM-Newton observations 
during the flare F2 from a Gem. The temperature and smission 
measure are corresponding to the hottest component of the best 
fit 3T-plasma model. 



were found to vary during the flare. The Z reached a value 
of 0.22 ± 0.01.^0 at flare peak and thereafter decreased to a 
value of O.16±O.O1Z0 at the bottom of the flare light curve. 
In the flare F7, temperature was peaked before the emis- 
sion measure. During the flare F7 temperature and emission 
measure were increased by a factor of ~ 3 and ~ 5.5, re- 
spectively. The X-ray luminosity reached a value of lO'^^'^ 
erg s~^, which is ~ 2 times more than that of the quiescent 
state. 



4.5 Loop modeling 

In order to trace the coronal structure responsible for the 
flares, we employed time dependent hydrodynamic model of 
Reale et al. (1997) for a single loop and a simplified approach 
as given by Haisch (1983). The details of hydrodynamic flare 
model (Reale et al. 1997, 2004 and Reale 2007) are given in 
Pandey & Singh (2008). This model includes both plasma 
cooling and the effect of heating during flare decay. The 
empirical formula for the estimation of the unresolved flaring 
loop length has been derived as 



4.4.5 Flare F7 from EI En 

Fig. [12] shows the spectral evolution of the flare F7 of the 
star EI Eri. Similar to flares observed from V711 Tan, we 
have flt the 4T plasma model keeping flxed 3T model pa- 
rameters from quiescent state. Nh and Z were left free to 
vary during the spectral fltting. We did not observe any sig- 
nificant variation in Nh during the flare F7. As shown in 
Fig. [12] and table [6] Z, temperature and emission measure 



Td-s/Err. 



3.7 X 10-4/(C) 



/(C) > 1 



(1) 



where Lhmd is loop half-length in unit of cm, Emax is the 
maximum temperature in unit of K, r^j is the decay time 
of the light curve, and f{C,) is a non-dimensional correction 
factor larger than one. C, is the slope of the decay path in the 
density-temperature (n-T) diagram (Sylwester et al. 1993) 
and is maximum (~ 2) if heating is negligible and minimum 
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Table 6. Spectral parameters for each temporal segment of the 
flares as marked by vertical lines in Fig. [T] 



Object (LS) FS Nh 



kT 



Lx xi (DOF) 



UZ Lib (Fl) Dl 
D2 
D3 
D4 
D5 



A And(F3)' 



D9 

{F2)* Ri i.elJJ-j 

R2 1.7+Jol 

R3 2.01°;^ 

Dl 2.2lg;2 

r)2 2 7+^-^ 
^-'-0.2 

D3 2.6l°-2 

D4 2Ati-l 

D5 2Atl-l 

D6 2.4lg;^ 

D7 2.2lg;2 

D8 2.2^0-2 

D9 2.4tg;^ 

Dio 2.3j:°;2 

Dll 2.6j:°'2 
Rl 2.9lg;? 
R2 3.6^2 ? 
R3 2.9+°-^ 
R4 3.2+_i% 
R5 3.lt°-? 
3 + 0-6 



R6 2.9 



V711 Tau {F4) Rl 
Dl 

{F5) Rl 
R2 
Dl 
D2 
D3 
D4 
D5 
D6 

{F6) PF 
Rl 
R2 
R3 
Dl 
D2 
D3 
D4 
D5 
D6 
D7 

EI Eri (FT) Rl 
R2 
Dl 
D2 
D3 
D4 
D5 
AF 



+0.6 
■ 0.4 
. + 0.5 
' 0.5 
,+0.6 
0.6 
,+0.5 
■0.5 
, -,+0.5 
^^■^-0.5 
1 1 + °-^ 
.+0.4 
■0.4 
. + 0.3 
0.3 
+ 0.1 
0.1 
+0.2 
0.2 
+0.2 
■0.2 
1 2+''-2 
-, ,+0.2 
^■■^-0.2 
+0.2 
■0.2 
,+0.2 
■0.2 
,+0.2 
0.2 

+0.2 
0.2 
1 1 + 0-2 
, + 0.4 
0.4 
,+0.5 
■0.5 



LI"- 
1.6l 
2.oj 
2.0+ 



1-6: 
1.6+ 
1.1-1 
l.t-l 
l.t-l 



0.9+ 



LI"- 



0.8 

1.3 

< 0.8 

5+'^-* 
"■ = -0.4 

7+'^-3 
"■'-0.5 

9+'^-'' 
"■''-0.4 

-^■^-0.4 

6+'^-2 
"■°-0.2 



32+°-"^ 
"■■^■^-0.02 

29+° °^ 
^■•^^-0.02 
„ ,0 + 0.02 
"■^°-0.02 

"■•^"-0.02 
0.28^ 



, + 0.02 
■0.02 
,+0^02 
0^01 
27+°"^ 



0.27^ 



0.25 
0.26 
0.26 



0.02 
+ 0^02 
0^01 
+ 0^02 
0^02 
+ 0.02 
0.02 



0.24^ 



+ 0.02 
■0.02 

"■^■^-0.02 

"■^■^-0.02 

"■^^-0.02 

o.ir 



-0-6 



+ 0.01 
■0.01 
,+0.01 
■0.01 
+ 0.01 
■0.01 
,+0.01 
■0.01 
+ 0^02 
0^01 

"■^-'-o^oi 



0.091 
O.llj 

o.ioj 

0.11^ 



0.2123 
O.173I 
0.169j 
0.170+ 



0.191: 
0.188j 
0.184j 
0.174j 
0.170j 
0.165j 
0.1781 
0.217j 
0.216j 
O.215I 
0.182j 
0.174j 
0.162j 
0.165+ 



4 0+0- = 
^■"-0.4 
4 4+0-6 

4 9+0-8 

^■■'-0-6 

■^■"-0.5 

3 9+0-7 
■^■^-0.6 

4.6+ 



+ 1.2 
0.8 
+ 0.8 
0.6 
+ 1.7 
1.0 
+ 0.2 
0.2 
r, Q+0.2 
0.2 
r, 7+0-2 
0.2 
+ 0.1 
0.1 
+ 0.1 
0-1 
+ 0.1 
0-1 
+ 0.1 
0.1 
+ 0.1 
0.1 
+0.11 
0.08 
+0.08 
0.08 
+0.08 
0.08 
+0.09 
0.08 
+0.09 
0.09 
+ 0.08 
-0.08 
+0.6 
0.4 
+0.8 
0.7 
+0.5 
0.4 
„ 0+0.5 
•^■8-0.4 
+0.4 
0.4 
+ 0.2 
0.2 



3.73 
3.6j 
6.4j 
,.9+ 
■■71 
5.0j 
4.9I 
4.6l 
4.71 
4.4I 
4.37 
4.19 
4.00 
3.95 
3.86 
3.6 
3.4 
4.6 
3.7 
3.8 
3.4! 
2.9+ 



20h+°"^^ 
"■^"O-O.oil 
„ ,,0+0.012 
"■^^■^-0.012 

"■^^'-0.014 
o ,,-+0.013 
"■^^'-0.013 

234+°^°'-2 
"■^•^*-0.012 
n o,q+0.015 
"■^^"-0.010 
0.222^' 



+ 0.009 
0^009 
+ 0^008 
O^OOS 
+ 0.003 
0.004 
+ 0^005 
0.005 
+ 0^006 
0^006 
n ,«9 + 0.006 

n 184+0-004 

+ 0^006 
0^006 
+ 0^006 

o^ooe 

+ 0.007 
0.007 
+ 0.007 
0.007 
+ 0.010 
0.005 
+ 0.009 
0.010 
+ 0.011 
0.011 
+ 0.014 
0.014 
+ 0.014 
0.014 
+ 0.011 
0.011 
+ 0.010 
0.010 
+ 0.009 
0.009 
+ 0.011 
0.011 
+ 0.006 
0.005 



< 3.3 
3+1.4 
-0-9 
3+1.3 
-0.9 
-,+0.8 
-0.5 
5+0.6 
-0.3 
3+0-3 
-0.5 
,+0.3 
-0.4 
,+0.5 
-0.2 
3+0.3 
-0.2 
,+0.4 
-0.3 
5+0.2 
'-0.1 
,+0.4 
-0.2 
1+0.2 
-0.2 
, ^+0.2 
^■°-0.2 
1+0.2 
-0.2 
3+0.2 
-0.2 
3+0.3 
-0.2 
7+0.3 
-0.3 

^■°-0.2 

1 5+°^2 
^■°-0.1 
1 0+0.2 
^■•^-0.2 

^■^-1.1 
7+0.6 
-0.6 
,+0.7 
-0.5 
,+0.8 
-0.5 
3 o+^^i 

•^■"-0.8 
,+2.5 
-1.1 
,+0.9 
-0.4 
,+0.7 



3.93 

2.9I 

4.0j 

2.6j 

2.! 

2.33 

2.lj 

2.0j 

2.1^ 

1.8-t 

2.lj 

2.4+ 



2.43 
2.0j 
2.0j 
1.7+ 



3.7Z 
3.3j 
3.2+ 



3.1: 

1.6:: 

3.2"! 



81^3 



46l^ 
351^3 
33l^ 
3ll« 
3461^ 

7+7 



3571 
377 



3701° 
365^^ 
339^^ 
332j;* 
325t^ 
320l^ 
315±^ 
307j;^ 
292+5 
286^^ 
a+1 



a+1 



+ 2 
-1 
7+2 



6.1: 
9.5+ 



14 

17_^ 
, ,+0*7 
^■^-0.7 
+ 0.9 
0-9 
+3.2 
2.1 
27.7t^'^ 
36.7ii;^3 

20.9t|i^ 
16.7 



0.5 



.+2.1 
-2.1 
+ 1.S 
1.8 
,+ 1.0 
■0.9 
+ 1.4 
1.7 
, + 1.7 
1.6 

'+;"7 

25.it;;| 

,+ l.S 
1.7 
,+ 1.7 
1.7 
. + 1.9 
'-1.9 
, + 1.9 

■ 1.9 
, + 1.7 
'-1.7 

+ 2.5 
■2^3 
,+2.8 
■2.8 
.+5.6 
■3.7 
,+4.0 
■4.0 

■ + 2.S 
' 2.8 
,+3.6 

■3.6 
o ,+2.8 
'■^-2.8 
9.2t-° 
, + 1.5 

■ 1^4 



12.53 

11. oj 
12. ij 
19. 2j 
22.9+ 



20.93 
14.71 
11. ( 
11.33 

9.8j: 

10. 4j 
21. oj 
45. 7I 
31. Ol 
16.61 
12.9+ 



10.8^ 



17.98 
17.24 
16.20 
15.73 
14.28 
14.12 
13.13 
12.54 
12.29 
64.48 
65.40 
65.18 
64.12 
62.80 
61.09 
58.96 
56.00 
54.05 
52.49 
51.34 
49.86 
48.11 
46.38 
3.74 
4.15 
4.47 
4.58 
4.62 
4.61 
6.1 
6.5 
7.84 
10.52 
10.75 
10.67 
9.79 
9.11 
8.39 
7.88 
7.29 
7.44 
8.28 
8.80 
8.96 
8.51 
8.03 
7.55 
7.28 
7.14 
7.02 
10.37 
13.64 
11.89 
10.09 
8.85 
8.25 
7.97 
8.55 



0.99(315) 
1.09(311) 
0.97(289) 
1.07(291) 
1.05(283) 
0.95(278) 
0.91(267) 
1.08(260) 
1.09(89) 
1.07(1428) 
1.05(1345) 
1.17(1332) 
1.11(1367) 
1.15(1385) 
1.18(1361) 
1.10(1354) 
1.13(1308) 
1.07(1288) 
1.17(1263) 
1.07(1249) 
1.16(1223) 
1.05(1199) 
1.16(1166) 
1.22(513) 
1.15(551) 
1.06(576) 
1.13(587) 
1.07(582) 
1.20(587) 
0.95(511) 
1.12(523) 
1.32(246) 
1.27(285) 
1.14(250) 
1.35(290) 
1.10(282) 
1.22(290) 
1.35(315) 
1.36(338) 
1.20(627) 
1.11(479) 
1.10(462) 
1.11(495) 
1.11(527) 
0.93(476) 
1.13(457) 
1.09(440) 
1.26(430) 
1.16(424) 
0.96(423) 
0.98(322) 
1.02(319) 
1.02(294) 
0.90(323) 
1.14(341) 
1.12(312) 
0.93(352) 
1.16(477) 



Note: LS stands for light curve segments, FS stands for "Flare Seg- 
ments", iVjf is in unit of 10*^0 cm-"^, abundances Z is in solar unit, 
temperatures (kT) are in keV, emission measures (EM) are in 10^*^ cm""^, 
and X-ray luminosity (Lx) is in 10-^0 erg s~^. xi^ is the minimum reduced 

and DOF stands for degrees of freedom. 
* Spectral parameters of hot component in APEC 3T fit 




5 10 15 20 25 

Time (in ks since ,ID 2451936.^406) 



Figure 9. Similar to Fig. [8] but for flare F3 from star A And 




10 15 20 25 30 

Time (in ks since -ID 2451963.0728) 

Figure 10. Similar to Fig. |8] but for the flare F5 from star V711 
Tau 



(~ 0.5) if heating dominates the decay. The correction fac- 
tor fiC,) and loop maximum temperature are calibrated for 
the XMM-Newton EPIC spectral response and are given as 
(Reale 2007): 

0.51 



/(C) 



+ 1.36 



(for 0.35 < C < 1.6); (2) 



C- 0.35 

The average temperatures of the loops, usually lower than 
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Figure 11. Similar to Fig. [8] but for tlic flare F6 from star V711 
Tau 
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Figure 12. Similar to Fig.[8]but for the flare F7 from star EI Eri 



the real loop maximum temperatures, are found from the 
spectral analysis of the data. The loop maximum tempera- 
ture for EPIC instruments is calibrated as 



(3) 



The loop length derived from this method is purely 
based on the decay phase of the flare. Alternatively, Reale 



(2007) derive the semiloop length from the rise phase and 
peak phase of the flare as; 



Lhn 



rpbl-i 

9.5 X 10 tM^p^ 



(4) 



where L-jimr is the loop half-length in unit of cm, Tmai is 
the maximum temperature in unit of K, Tm is the tempera- 
ture at density maximum and t m is the time in unit of s at 
which density maximum occurs. The time of the maximum 
emission measure is a good proxy for Tm - Haisch (1983) sug- 
gested a simplified approach which is based on quasi static 
radiative and conductive cooling during the earliest phases 
of flare decay. Given an estimate of two measured quantities, 
the emission measure {EMave) and the decay time scale of 
the flare {Td ) Haisch's approach (Haisch 1983) leads to the 
following expression for loop length (Lha) 



Lha = 5x 10" 



(5) 



Once the length of the flaring loop has been derived, 
it is possible to infer additional physical parameters of the 
flaring plasma. These parameters could be electron density 
(jie), flaring volume {V), maximum pressure (p), heating 
rate per unit volume {Eh) and minimum magnetic field (B) 
to confine the plasma. For details, we refer to our previ- 
ous paper (Pandey & Singh 2008) . The results of the model 
parameters are given in Table [7] and explained below. 



4.5.1 UZ Lib: Flare Fl 

The maximum temperature was determined using the equa- 
tion [3] and found to be 128 ± 21 MK for the flare Fl. Fig. 
1131 (a) shows the evolution of the flares in n-T diagram with 
a least-square fit to the decay phase, where the EM'^'^^ has 
been used as a proxy of density. The path of n-T diagram is 
shown by dotted line with arrows. The continuous straight 
line shows the best linear fit to the corresponding data dur- 
ing the decay phase, providing a slope The value of the 
indicates the presence of sustained heating during the de- 
cay phase of the flare Fl. The value of ^ for the flare Fl 
is less than the lower asymptotic value for which the equa- 
tion [T] can be applied. Therefore, an upper limit for the loop 
length can be derived by using the value of ^ at the high 
extreme of the error bar. The resulting loop length using 
the hydrodynamic loop model based on decay phase are es- 
timated to be < 18.2 X 10^° cm (~ 0.2R*) for the flare Fl. 
The loop lengths derived from hydrodynamic method based 
on decay phase is found to be ~ 1.3 times smaller than 
that derived from Haisch approach (see Table 7). The esti- 
mated loop lengths were found to be much smaller than the 
pressure scale heigh10 [hp ~ 5 x 10^^ cm). The maximum 
pressure in the loop at the flare peak, estimated from the 
loop scaling law (Rosner, Tucker & Vaiana 1978) was found 
to be > 4300 dyne cm"'^ for the flare Fl. Assuming that 
the hydrogen plasma is totally ionised (p — 2nekTmax) , the 
maximum plasma density in loop at the flare peak is esti- 
mated to be > 1.2 X 10^^ cm~^ for the flare Fl. Using the 
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26.70 26.80 26.90 27.22 27.24 27.26 27.28 27.30 

0.5 log EM (cm"^) 0.5 log EM (cm"^) 



(a) Flare Fl (b) Flare F2 





26.4 26.6 26.8 

0.5 log EM (cm"^) 

(e) Flare F7 



Figure 13. The density-temperature (n-T) diagram during the decay phase of the flares, where EM^/^ has been used as a proxy 
of density. Symbols for the each flares are given at the bottom right corner. Straight lines represent the best linear fit to the data 
corresponding to the decay phase. C, is slope of n-T diagram. The dotted lines with arrows show the path of n-T diagram. 



observed peak emission measure the loop volume was com- 
puted to be of the order of < 5.4 x 10^^ cm^. A hint for the 
heat pulse intensity comes from the flare maximum temper- 
ature. By applying the loop scaling laws (Rosner, Tucker 
& Vaiana 1978; Kuin & Martens 1982) and loop maximum 
temperature, the heating rate per unit volume {Eh) for the 
flare Fl would be > 0.73 erg cm^'^ s^^, respectively. From 
the derived pressure of the flare plasma, the minimum mag- 
netic field required to confine the plasma should be 328 G. 



4.5.2 (J Gem: Flare F2 

The maximum temperature of the flare F2 was derived to 
be 177 ± 6 MK, which is the highest among all the observed 
flares. The slope of n-T diagram indicates that the sustained 
heating was negligible during the decay phase of the flare. As 
seen from Fig. [T] (c) and Fig. |8] the quiescent states before 
and after the flare were not observed. Therefore, the larger 
value of the C, could be due to incomplete flare evolution. The 
value of ^ for the flare F2 is outside the domain of the valid- 
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Table 7. Parameters derived for flares. 



Object(FN ^) (1) 
Parameters (^) 


TTr7T "1 / "n -1 \ 

UZLib(Fl) 


CT Gem(F2) 


A And(F3) 


V711 Tau (F4) 


V711 Tau(F5) 


V711Tau(F6) 


El Eri(F7) 


c 


(2) 


0.21 ±0.38 


2.17 ±0.33 






0.65 ±0.33 


0.97 ±0.41 


0.6 ±0.5 


^max 


(3) 


128 ± 21 


177 ±6 


121 ± 24 


98 ± 32 


101 ± 17 


59 ±5 


129 ± 30 


Tm 


(4) 




133 ±4 


70 ±7 




61 ± 11 


56 ±5 


93± 17 


EMave 


(5) 


50.1 ±3.6 


336 ±6 


11.8 ± 1.7 


4.2 ± 0.8 


23.8 ± 1.0 


14.8 ±0.6 


20.7 ±1.3 


tM 


(6) 




> 11250 


25850 ± 2350 




5900 ± 900 


5750 ± 750 


1900 ± 500 




(7) 


< 18.2 


< 64.3 






< 10.6 


6.5 ± 1.7 


< 5.5 


^hmr 


(8) 




> 25.5 ±4.1 


79.7 ± 18.3 




11.9 ±3.2 


4.6 ±0.8 


3.9 ± 1.5 


Lha 


(9) 


23.4 ±0.7 


51.8 ±0.4 




2.8 ±0.3 


9.9 ±0.2 


7.5 ±0.1 


5.2 ±0.1 


P 


(10) 


> 4.3 


> 3.1 


0.8 ±0.5 


13± 12 


> 3.6 


1.2 ±0.4 


> 14 


< 


(11) 


> 1.2 


> 0.6 


0.24 ± 0.07 


4.6 ± 1.6 


> 1.3 


0.7 ±0.2 


> 3.9 


V* 


(12) 


< 54 


< 919 


293 


0.3 


< 22.3 


49.8 


< 2.9 


Eh 


(13) 


> 0.73 


0.18 


0.03 


12.3 


0.94 


0.38 


7.8 


B 


(14) 


328 


279 


141 


571 


300 


173 


593 



(1) FN is flare name, (2) f is slope of decay path in density-temperature diagram, (3) Tmax is the maximum temperature in the loop at 
the flare peak in unit of lO'' K based on spectral flt and equation l(3|, (3) Tm is the temperature in the loop at the density peak in unit 

of lO'^ K based on spectral fit and equation JSj, (4) EMave. is the average emission measure during the flare in unit of 10^^ cm~^, (5) 
Im is the temperature at which density peaks in unit of seconds (7,8,9) Lhmd, Lfimr and L^a are lengths of the flaring loop in unit of 

10^" cm and are determined using equations |[T}, l(4]|, and l(5j, (10) p is the maximum pressure in the loop at the flare peak in unit of 
10^ dyne cm and. derived from R,XV law, Tmax — 

1.4 X 103(pL)^/^ (Rosner, Tucker & Vaiana 1978), (11) rxe is the maximum 
electron density in the loop at the flare peak in unit of 10^^ cm~^, assuming that hydrogen plasma is totally ionised (p = 2ne.kTmax) , 
(12) V is the volume of flaring plasma in unit of lO''" cm^ and estimated using the equation EM = n^V, where EM is the peak 
emission measure. (13) E^ is the heating rate per unit volume at the flare peak in the unit ergs s~^ cm~^ estimate from the loop 
scaling law. Eh ^ lO'^Tj^^^L-'^ . (14) B is the minimum magnetic field necessary for confinement in Gauss and can be simply 

estimated as p = B^ /8tt. 

(*) The pressure as derived here is a maximum value, appropriate for the equilibrium condition, which can be reached only for very 
long-lasting heating. The density and volume values derived thereby are therefore be treated as upper and lower limits, respectively. 



ity of method, therefore, loop length of < 64.3 x 10 cm was 
derived by using ^ — 1.6 (the maximum allowable value). Al- 
ternatively, the loop length can be derived using the equa- 
tion Q. As seen in Fig.[TJc), the entire rise phase of the flare 
F2 was not observed, therefore, the value of tM was deter- 
mined to be in a lower limit. Thus, the loop length derived to 
be > 25.5 ±4.1 x lO^'' cm, which is ~ 2.5 times smaller than 
that derived from the hydrodynamic decay model. The loop 
length derived from hydrodynamic method based on peak 
and rise phase is also found to be ~ 2 times smaller than 
that from Haisch approach (see Table [T]). The derived loop 
length from equation ([T)) may be over estimated as the value 
of C, is at the upper extreme end. However, loop length de- 
rived from hydrodynamic method based on rise phase may 
be under estimated as the rise phase of flare was not ob- 
served completely. Loop length derived from each method is 
less than the stellar radius (i.e. L ~ 0.3 - 0.8R*). For further 
parameter estimation, we have used the loop length derived 
from hydrodynamic method, which is much less than the hp 
(=7 X 10^'^ cm). The estimated loop parameters are given in 
Table El 



4.5.3 X And: Flare F3 

In the case of the flare F3, the loop length could be derived 
only from the hydrodynamic model based of the rise and 
peak phases (see equation|4| as no decay phase was observed. 
The derived loop length (=79.7 ± 18.3 x 10^° cm 1.5 R*) 
is less than the pressure scale height (=5 x 10^"^ cm). The 
maximum temperature was estimated to be 121 ± 24 MK. 
The maximum pressure, density and volume are estimated 



to be 800±500 dyne cm~^ 2.4±0.7x 10^" cm"^ and 2.93 x 
10'^^ cm'', respectively. The heating rate per unit volume 
was estimated to be smaller than that estimated for other 
observed flares. 



4.5.i V711 Tau: Flares F4, F5 and F6 

For the flare flares F4, F5 and F6 the loop maximum tem- 
peratures were found to be 98 ± 32, 101 ± 17 and 59 ± 5 
MK, respectively. The maximum temperature observed dur- 
ing the flare F6 is the lowest among all the other observed 
flares. The slopes of decay path of n-T diagram for the flares 
F5 and F6 (see Figs. [T3b and I13t i) indicate the presence of 
a strong sustained heating showing that the observed decay 
paths were driven by time evolution of heating process. For 
the flare F4, loop length could be derived only from Haisch 
approach and the derived loop length was 2.8 ± 0.3 x 10^'^ 
cm (~ 0.2/O.lR*). For the flare F5, the loop lengths derived 
from each of the methods are consistent within la level. 
Adopting the loop length derived from hydrodynamic decay 
method (i.e. <0.8/0.4 R*) and the peak emission measure of 
37^5 X 10^^ cm~^ for the flare F5, the maximum pressure, 
density and volume were estimated to be > 3600 dyne cm~^, 
> 1.3 X 10" cm"^ and 2.23 x 10^^ cm^, respectively A 300 
G of minimum magnetic field was estimated to confine the 
plasma in the flaring loop for the flare F5. For the flare F6, 
the loop length derived from hydrodynamic decay method 
is within ~ la level to that derived from hydrodynamic rise 
method. However, loop length derived from Haisch approach 
is more than that from hydrodynamic methods. Loop pa- 
rameters for the flare F6 were derived by adopting the loop 
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length of 6.5 ± 1.7 x 10^° cm (~ 0.24/0. 5R*) and are given 
in Table [T] The loop length derived for each flare is less then 
the pressure scale height (~ 4 - 



4.5.5 EI En: Flare F7 

Fig. I13f e') shows the n-T diagram for the flare F7, where the 
slope of 0.6 ± 0.5 implies the presence of sustained heating 
during the decay. The loop length based on equation [T] was 
derived by using the upper extreme value of ^ = 1.1 as ^ 
at lower extreme of the error is outside the validity of the 
method. Thus the loop length was derived to be < 5.5 x 10^° 
cm (~ 0.23R*), which is similar to that derived from equa- 
tion [l] The loop lengths derived from hydrodynamic meth- 
ods based on rise and peak phases is smaller than that de- 
rived from the Haisch approach (see Table [T]). The derived 
loop length is smaller than the pressure scale height. The 
maximum temperature of the loop was derived to be 129±30 
MK. Using the loop length of 5.5 x 10^" cm and peak emis- 
sion measure of 4.57 x 10^'' cm"'', the maximum pressure, 
density and volume of the flaring plasma were estimated to 
be > 1.4 X 10* dyne cm'^ > 3.9 x lO" cm"^ and 2.9 x 10^° 
cm^, respectively. A minimum magnetic field of 593 G was 
required to confine the plasma at the peak of flare. 



5 DISCUSSION AND CONCLUSIONS 

We have carried out a detailed analysis of seven X-ray flares 
observed with XMM-Newton from five cool subgiants and 
giants in RS CVn binaries. Except the flare F4, the decay 
time (> 1 h) and the total energy released (1.2 — 8.0 x 10^^ 
erg s~^) indicate that these flares belong to a class of long 
decay flares (Pallavicini et al. 1990). However, excluding the 
flare F2 of a Gem, a strong sustained heating was found to 
be present during the decay of all flares. The strong sus- 
tained heating during the flare decay indicates that the de- 
cay light curves is dominated by the temporal heating proflle 
and not by free decay of the heated loop. Therefore, these 
flare can not be classifled as long decay flares but may be 
classifled as an arcade. The decay of these observed flares 
are smaller than that of few pre-main-sequence stars, where 
the decay time was found upto few days (Favata et al. 2005, 
Getman et al. 2008). But, the observed flare decay time is 
comparable or more than that of the flare observed in the 
late type G-K dwarfs, other RS CVns and the flare dMe 
stars (e.g. Pandey & Singh 2008, Schmitt 1994, Osten & 
Brown 1999). Fig. [14] shows that correlation between the Td 
and Tr for the flares from RS CVns and G-K dwarfs. The 
straight lines show the linear regression flt in the form of 
Td oc T,"-6^±0-^« for RS CVns and n oc r° «*±«-20 for G-K 
dwarfs. These relations are similar to the relation predicted 
for standard magnetic reconnection model i.e. Td oc r°'^ ( 
Petschek 1964, Shibata & Yokoyama 1999,2002, Imanishi et 
al. 2003), where it is assumed that the decay phase is dom- 
inated by radiative cooling. 

There were two flares, F2 from a Gem and F5 from 
V711 Tau that showed a sustained level of approximately 
constant emission at the peak of the flare. In the past such 
scenario was seen only for a small number of flares from 
V711 Tau, o-^ CrB and AT Mic (e.g. Agrawal et al. 1986, 




iog(V 



Figure 14. Plot between log(r,.) and log(Ttj) for flares from a 
sample of RS CVns (solid circles) and G-K dwarfs (solid squares). 
The lines represents the best fit log-linear correlations. The values 
of Td and Tr for G-K dwarfs are taken from Pandey& Singh (2008). 



Osten & Brown 1999, Raassen et al., 2003). The appear- 
ance of plateaus could be due to the reconnection process, 
where reconnection balances the loop expansion to sustain 
a constant luminosity. 

During the flares observed in a sample of RS CVns, 
the peak luminosities were in the range lO'"''^ < Lx < 
10^^'* erg s"'^ indicating that these flares are more pow- 
erful than the flares observed from dwarf stars ( lO'^*'^ < 
Lx< 10^°-^; Pandey & Singh 2008) but are less powerful 
than the flares from PMS stars (10^^ < Lx< 10^^ erg s"^; 
Getman et al. 2008). This conflrms the previous finding that 
the large fiares on active RS CVn binaries involve small en- 
ergies in comparison to the total stellar luminosity. However, 
few fiares have been reported earlier from dMe stars, where 
peak X-ray flare luminosity reached upto 30 per cent of bolo- 
metric luminosity (e.g. GINGA EQ 1839.6-1-8002 -Pan et al. 
1997; EV Lac -Favata et al. 2000c; AB Dor - Maggio et al. 
2000) . The total energy released during the X-ray flares ob- 
served in our sample was found to be of the order of IQ'^^"''^ 
erg. This implies that these flares are lO'^"* more energetic 
than the flares observed in dMe and G-K dwarfs (Pandey & 
Singh 2008, Pallavicini 1990), but are as energetic as flares 
observed in some of pre-main-sequence stars (Favata et al 
2005; Getman et al. 2008) . The energy released during flares 
are also comparable to that from the superflare on the RS 
CVn binary H Peg (10^** **; Osten et al. 2007). 

It appears that the quiescent coronae of the stars UZ 
Lib, (T Gem and A And are represented by two temperature 
plasma. However, the quiescent coronae of V711 Tau and 
EI Eri are represented by three temperature plasma. The 
coolest and the hottest temperature of quiescent coronae 
were found to be ~ 0.3 - 0.8 and ~ 1.2-2.5 keV, respec- 
tively. These temperatures are more than that derived for 
RS CVns by Dempesey et al. (1993) using ROSAT PSPC 
data. They found the bimodal distribution of temperatures 
centered near 0.17 and 1.38 keV. The quiescent coronae of 
RS CVn -type stars appears to be hotter than the quiescent 
coronae of cool dwarf stars. The time resolved spectroscopy 
of these flares shows the coronal spectrum during the flare 
can be represented with a IT model, with quiescent state 
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taken into account as a frozen background contribution. In 
the case of a Gem, the Nh was found to vary along the rise 
and decay phases of fiare hght curve. Enhancement of hydro- 
gen column density, associated with flaring events, have been 
observed in the past in Proxima Cen (Haisch et al. 1983), 
V773 Tau (Tsuboi et al. 1998), Algol (Ottmann & Schmitt 
1996; Favata & Schmitt 1999), and LQ Hya (Covino et al. 
2001), and usually interpreted as due to Coronal Mass Ejec- 
tion (CME). As shown in Fig |5] the hydrogen column density 
increased as flare rise, peaked and then showed decreasing 
trend as flare decayed. Such a behavior is compatible with 
the picture of a cloud of material being ejected from the flar- 
ing region at the beginning of the event, and subsequently 
thinning out as it expands and cools into the space. So- 
lar observations indicate a tight relationship between CME 
and flares (Vrsnak et al. 2004, Goplaswamy 2010). How- 
ever, there are large number of solar flares without CME 
associations (Gopalswamy et al. 2009). Rodono et al. (1999) 
suggesting that extra-absorption might occur in some stars 
due to the presence of neutral hydrogen in circumstellar en- 
vironment. 

The metal abundances were also found to vary along 
the flare light curves in the flares observed from a Gem, 
V711 Tau and EI Eri (see Figs. [Hl[Tni [Hand [HI). Using the 
RGS data of flares F2 from a Gem, Nordon & Beher (2008) 
show that abundances enriched during the flare. A change 
of the metal abundance during a strong flare was detected 
for the RS CVn like stars II Peg (Mewe et al. 1997), UX Ari 
(Gudel et al. 1999) and cr^ Crb (Osten et al. 2003), Algol 
(Ottmann & Schmitt 1996; Favata & Schimtt 1999), dwarf 
stars EV Lac (Favata et al. 1998; 2000a), CN Leo (Liefke et 
al. 2010) and AB Dor (Giidel et al. 2001), and premainse- 
quence star V773 Tau (Tsuboi et al. 1998). In these cases 
abundances vary with an initial rise from a quiescent state 
value of Z ~ 0.3^0 to Z = IZq, and a subsequent decay to 
the pre-flare value. Our results are in general consistent with 
the previously observed behavior. In case of a Gem the metal 
abundances go up as material is ejected, as is seen in the si- 
multaneously increasing Nh- The most likely interpretation 
for the abundances variations observed during the flares is 
that at the beginning of the event fresh chromospheric ma- 
terial is evaporated in the flaring loops, thus enhancing its 
abundances; once the material is brought in the flaring coro- 
nal structure the fractionation mechanism which is responsi- 
ble for lower abundances observed in coronal plasma would 
start operating, bringing the coronal abundances back to 
quiescent state value (see also Favata & Micela 2003). 

For most of the flares, the temperature and the emis- 
sion measure show well-defined trends i.e. the changes in the 
temperature and the EM are correlated with the variations 
observed in the light curves during the fiares. Similar trends 
were also seen during the fiares detected in the cool pre- 
main-sequence and main sequence stars (Favata et al. 2005; 
Giardino et al. 2006; Pandey & Singh 2008). The peak flare 
temperature was found in the range of 59 - 177 MK for all the 
flares. These values are similar to those observed from ac- 
tive stars like pre-main-sequence stars {Tmax ~ 80-270 MK; 
Favata, Micela & Reale 2001; Favata et al. 2005; Getman et 
al. 2008), Algol (T„,ax ~ 70-150 MK; van den Oord & Mewe 
1989, Favata et al. 2000c), AB Dor (Tmax ~ 140-170 MK; 
Maggio et al. 2000), other RS CVns (T,„„^ ~ 25-120MK, 
Giidel 2004) but are more than those found in the flares 



from the dMe star and G-K dwarfs (~ 10-50 MK; Reale & 
Micela 1998; Favata et al. 2000b; Briggs & Pye 2003; Pandey 
& Singh 2008 ). The flare peak temperatures in our sample 
were also found to be comparable or more than that de- 
rived for a superflare on II Peg (118-152 MK; Osten et al. 
2007). In these flares either both the emission measure and 
the temperature peaked simultaneously or the temperature 
evolved before the emission measure. A similar delay is of- 
ten observed both in the solar and stellar flares (Sylwester 
et al. 1993; van den Oord & Mewe 1989; Favata et al. 2000a; 
Maggio et al. 2000; Stelzer et al. 2002). A delay between 
the temperature and the emission measure peak indicates 
a coherent plasma evolution and therefore a flare occurring 
inside a single loop, or at least the presence of a dominant 
loop early in the flare. 

For all the flares, the estimated maximum electron den- 
sity under assumption of a totally ionized hydrogen plasma 
was found in the order of 10^°~^^ cm^'' . This is comparable 
to the values expected from a plasma in coronal conditions 
(Landini et al. 1986). 

The observed flares are modeled using the hydrody- 
namic model based on the decay time and the rise time, 
and Haisch's approach. For the flares F5, F6 and F7, the 
loop lengths derived from the hydrodynamic model based 
on the decay time are found to be consistent within la to 
that from the rise time and peak phase method. In the case 
of flare F2 of a Gem, however, the loop lengths from the two 
methods were found to be inconsistent. This could be due 
to the involvement of several different coronal loops during 
the decay of the flare or the heat pulse triggering the flare 
may not a top hat function as is assumed normally. During 
the flares Fl, F5, F6 and F7, the decay paths are driven by 
sustained heating, which indicate that a flare is involving 
a loop arcade or is a two-ribbon flare. It may then occur 
that the flare begins in a loop and propagates to other loops 
only progressively, i.e. at later times. Therefore, we may ob- 
tain a signiflcant discrepancy between the estimation made 
from the rise phase and the one from the decay phase. Gen- 
erally, the loop lengths derived from Haisch approach were 
found be more than that derived from hydrodynamic mod- 
els. However, in most of the cases difference between the 
loop lengths obtained from three methods is well within la 
level. In the case of flare Fl, the loop length derived from 
Haisch approach are signiflcantly larger than that derived 
from the hydrodynamic model. Covino et al. (2001) have 
applied one version of the Haisch's approach to a number of 
X-ray flares on the stars. They report that the loop lengths 
of these flares derived from the simplified Haisch approach 
is in general agreement with the method based on detailed 
hydrodynamic calculations. However, Favata et al. (2001) 
showed that the large loop sizes obtained by quasi-static 
analysis are in reality likely to be significantly larger than 
the hydrodynamic model, typically by factors of > 3. 

The derived loop lengths presented in Table 7 can be 
compared to the loop lengths derived from other stars. The 
present analysis shows that RS CVn-type binaries have 
larger coronal structure which are comparable to or less than 
their radii. The log of loop length (log L) in our sample 
ranges from 10.4 to 11.9 . These loop lengths are similar 
to or more than that observed for the fiares from other ac- 
tive RS CVns like binaries e.g. UX Ari (log L = 11.0-11.7 ; 
Franciosini et al. 2001, Giidel et al 1999), a^ Crb (log L = 10- 
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10.17 ; Osten et al. 2003, van den Oord et al. 1988, Agrawal 
et al. 1986), II Peg (log L=10.46-ll ; Mewe et al. 1997, Doyle 
et al. 1991, Tagliaferri et al. 1991). However, these flares are 
comparable to the flares observed from Algol (log L = 10.3- 
11.9 e.g. Favata et al. 2000c, van den Oord & Mewe 1989, 
White et al. 1986). If we compare these loop lengths with 
the flares from much younger stars, it was noticed that the 
flaring structures are larger in RS CVn binaries than that 
from the main-sequence stars (log L =9.8-10.6; Pandey & 
Singh 2008), and are comparable or smaller than the flaring 
structures from pre-mainsequence stars (log L~10-12; Fa- 
vata et al. 2005, Getman et al. 2008). This imphes that RS 
CVns have large coronal magnetic activity. The high activ- 
ity level of RS CVns is generally due to their fast rotation, 
itself a consequence of tidal locking of the stellar rotation 
period with the orbital period of binary. 

The large sizes of flaring loops (larger than the com- 
ponent star) indicate that the coronal structures in active 
binaries are large, and would typically span both stars. For 
active RS CVn binaries, Uchida & Sakurai (1983) suggested 
the idea of inter-binary loops, with magnetic structures link- 
ing the two stars. We have compared the loop length with 
the distance between binary components and found that the 
flaring loop lengths are much smaller then the binary sepa- 
ration for the star V711 Tau (~ 3R* for circular orbit; R*is 
radius of more active component in binary system; Fekel 
1983). For the largest flare from V711 Tau (F6), the loop 
length is ~ 13 per cent of binary separation. In the case 
of flare F7 from EI Eri, the loop length is one tenth of bi- 
nary separation (~ 3 R* for circular orbit; Washuettl et al. 
2009). For other binary stars, we could not compare the loop 
lengths with binary separations as the separations between 
two components are not known, only the semi-major axis of 
one component is known (see Eker et al. 2008). The origin 
of flares in the close binary systems having starspots are in- 
terpreted in terms of the reconnection of the magnetic flux 
tubes of one component with the emerging and submerging 
pairs of spots on the other component (Uchida & Sakurai 
1985). However, the flares from V711 Tau and EI Eri appear 
to originate very close to the stars and extended to only a 
small fraction of binary separation. In the literature, it has 
been also noticed that corona of even the most active stars 
is quite compact and does not extend to large (larger than 
stellar radius) distance from the star itself (e.g. Brickhouse 
et al. 2001). 

At the flare peak, the maximum X-ray luminosity must 
be lower than the total energy rate [H = Eh-V) to sat- 
isfy the energy balance relation for the flaring plasma as a 
whole. The rest of the input energy is used for thermal con- 
duction, kinetic energy and radiation at lower frequencies. 
For the flares under study the maximum X-ray luminosity 
observed is about 40-60 per cent of H. These values are more 
than those reported for the solar flares, where the soft X-ray 
radiation only accounts for less than 20 per cent of total en- 
ergy (Wu et al. 1986). In comparison, the fraction of X-ray 
radiation to the total energy has been found to be 15 to 35 
per cent for the G, K and M dwarfs (Favata et al. 2000a,b; 
Reale et al. 2004, Pandey & Singh 2008). 

Applying loop scaling law V = 27r/3^L^, and if the de- 
tected flares are produced by single loop, their aspect ratio 
(/3) were estimated in the range 0.1 to 0.2 for all the flares. 
Similar aspect ratio was also observed for solar and other 



stellar coronal loops for which typical values of /3 are in the 
range of 0.1 - 0.3 (e.g. Shimizu 1995). 

We have derived the flare rates for RS CVn-types and 
G-K dwarfs. Following Stelzer et al. (2000), the flare rate 
(F) is estimated as 



F ■ 



t7 Td 



(6) 



where r = r J + r!^ and = ^ u^^ + are the mean 
value and uncertainty in e-folding decay and rise times de- 
rived with ASURV's Kaplan-Meier estimator, N is number 
of flares, Td is total duration of observations. For a sample 
of RS CVn-type stars, the fr and fd are determined to be 
8709 ± 1341 and 4168 ± 917, respectively. Total number of 
flares observed from the sample of RS CVn are 7 and total 
duration of the observations is 462 ks. Thus the flare rates 
for a sample of RS CVn stars was found to be 19±2 per cent. 
Similarly, using our previous paper (Pandey & Singh 2008) 
the flare rate for the G-K dwarfs was found to be 18.5 ± 0.3 
per cent. The flare rates for pre-mainsequence stars in p Oph 
star forming region and Taurus molecular cloud were deter- 
mined to be 15 ± 1 and 9 per cent, respectively (Imanishi 
et al. 2003; Stelzer et al. 2007). Surprisingly, we obtained 
the larger flare rates for RS CVns and G-K dwarfs. Gener- 
ally, it was shown that the younger stars tend to have more 
frequent flaring activity (see Imanishi et al. 2003). Larger 
flare rate for RS CVns with respect to young stars may be 
result of longer rise and decay times, and limited number of 
observations. 
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